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The  behavior  of  an  atomic  impurity  coupled  to  conduction  electrons  has  become 
one  of  the  paradigms  of  condensed  matter  physics.  Competition  between  on-site 
Coulomb  interaction  and  band  hybridization  produces  the  Kondo  effect:  a  crossover 
from  weak  to  strong  coupling  between  the  localized  and  band  electrons  below  the 
Kondo  temperature,  Tk-  The  study  of  the  Kondo  effect  has  been  limited,  however, 
by  the  nature  of  the  impurity  system.  Since  it  is  a  daunting  task  to  drive  the  host 
metal  out  of  equilibrium,  it  is  the  equilibrium  properties  of  Kondo  impurities  that 
have  been  explored.  ‘ 

In  this  paper  we  address  a  new  Kondo  system  in  which  non-equUihrium  is  rou¬ 
tinely  achieved,  namely  a  semiconductor  quantum  dot  weakly  coupled  to  its  leads. 
It  is  already  evident  that  Anderson’s  model*  for  a  Kondo  impurity  -  discrete,  in¬ 
teracting  levels  coupled  to  a  band  -  also  describes  quantum  dots.  Experimen¬ 
tally,  the  discrete  spectrum  of  a  single  dot  has  been  probed  by  transport’"*  and 
capacitance*  spectroscopy,  while  the  strong  on-site  Coulomb  interaction  is  observed 
in  Coulomb-blockade  conductance  oscillations.*'*’*  Theoretically,  Anderson’s  model 
has  provided  an  excellent  description  of  these  experiments  both  in  equilibrium,** 
and  non-equilibrium. However,  it  is  only  the  high  temperature  regime  that  has 
been  explored  experimentally,  while  it  is  below  Tk  that  the  Kondo  effect  emerges. 

Since  the  Anderson  Hamiltonian  describes  the  quantum  dot,  at  low  temperatures 


the  dot  must  behave  as  a  Kondo  impurity.  In  fact,  Glazman  &  Raikh"  and  Ng  k 
Lee^^  have  argued  that  at  zero-temperature  equilibrium  the  Kondo  resonance  in  the 
density  of  states  of  spin-degenerate  levels  will  always  lead  to  perfect  transparency 
of  the  quantum  dot  at  the  Fermi  energy.  In  contrast,  above  the  Kondo  tempera¬ 
ture,  resonant  tunneling  occurs  only  at  a  discrete  set  of  Fermi  energies.  Thus  the 
Kondo  effect  will  have  a  striking  experimental  signature  in  low-temperature  trans¬ 
port  through  a  quantum  dot.  Furthermore,  in  the  quantum  dot  system  the  leads 
coupled  to  the  dot  are  easily  biased  to  non-equilibrium  2md  the  dot  potential  can 
be  swept  continuously  with  a  gate.  Thus  new  physical  questions  which  were  not 
relevant  to  magnetic  impurities  can  be  raised.  In  particular,  what  happens  to  the 
Kondo  effect  out  of  equilibrium  Since  transport  measurements  on  single  quantum 
dots  require  significant  applied  bias,  this  question  is  of  immediate  importance. 

In  this  Letter  we  combine  several  approaches  (non-crossing  approximation,^* 
equations  of  motion, perturbation  theory,  variational  wavefunction  calculation^*) 
to  present  a  consistent  picture  of  low-temperature,  non-equilibrium  transport 
through  a  quantum  dot.  For  spin-degenerate  levels  at  equilibrium,  the  Kondo  pcak^^ 
in  the  density  of  states  at  the  chemical  potential  (Fig.  la)  leads  to  resonant  trans¬ 
mission  through  the  dot."'^^  A  voltage  bias  between  the  left  and  right  leads  causes 
the  Kondo  peak  to  split,  leaving  a  peak  in  the  density  of  states  at  the  chemical  po- 


tential  of  each  lead  (Fig.  lb).  We  find  that  the  amplitudes  of  the  split  Kondo  peaks 
are  suppressed  by  a  finite  non-equilibrium  lifetime.  This  lifetime  results  from  dissi¬ 
pative  transitions  in  which  electrons  are  transferred  from  the  high  chemical  potential 
lead  to  the  low  chemical  potential  one.  As  the  voltage  bias  is  increased,  the  lifetime 
decreases,  resulting  in  a  suppression  of  the  Kondo  peaks  and  thus  a  suppression  of 
the  differential  conductance  (Fig.  2a). Upon  application  of  a  magnetic  field,  the 
Kondo  peaks  shift  away  from  the  chemical  potentials  by  the  Zeeman  splitting,  but 
in  opposite  directions  for  each  spin  (Figs.  Ic  and  Id).  Interestingly,  therefore,  when 
the  chemical  potential  splitting  equals  the  Zeeman  splitting,  a  Kondo  peak  shifted 
away  from  one  chemical  potential  crosses  the  other  chemical  potential.  We  predict 
an  observable  peak  in  the  differential  conductance  at  this  crossing^*  (Fig.  2b). 

We  model  the  quantum  dot  and  its  leads  by  the  Anderson  Hamiltonian’ 

^  -I-  Yi  +  U'n,ni  +  53  (1) 

(T-.kfL.K  <'  9^,k^L.K 

where  cl^{Ckv}  creates  (destroys)  an  electron  with  momentum  k  and  spin  c  in  one  of 
the  two  leads,  and  c*{Cf)  creates  (destroys)  a  spin-cr  electron  on  the  quantum  dot. 
Since  we  are  interested  in  temperatures  smaller  than  the  orbital  level  spacing  in  the 
quantum  dot,  we  consider  only  a  single  pair  of  levels  on  the  dot  with  energies  6t  = 
eo+Ac/2andci  =  eo-Ae/2.  The  third  term  in  f  1)  describes  the  Coulomb  interaction 


between  the  two  localized  spins  which  we  take  to  forbid  double  occupancy (f/ 
oo),  while  the  fourth  term  describes  the  hopping  between  the  leads  aind  the  dot. 

Our  aim  is  to  calculate  the  current  through  the  dot,  7,  which  for  the  cawe 
of  proportionate  couplings  to  the  leads,  r‘(u;)  =  Ar"(a;),  where  = 

2ir£*€iKJ»)  expressed^®  in  terms  of  the  density  of  states, 

-iImG;(u/),  as 

-  />(“)!  f.M  [-;ImG;(u.)].  (2) 

In  Eq.  (2),  r„(w)  =  r‘(u;)r"(a/)/  [r‘(w)  +  r^(u;)],  and  G'{uf)  is  the  Fourier  trans¬ 
form  of  the  retarded  Green  function,  G^{t)  =  -tO(t)({c»(0>cJ(0)}). 

In  order  to  calculate  the  Green  function  G'(u/)  we  use  both  the  non-crossing 
approximation^^  and  an  equations-of-motion  method.*®'*  The  non-crossing  approx¬ 
imation  is  based  on  an  exact  mapping  of  the  infinite-!/  Anderson  Hamiltonian  (1) 
onto  a  slave-boson  Hamiltonian.  If  vertex  corrections  are  neglected,  the  propagators 
for  the  boson  and  the  fermion  degrees  of  freedom,  which  correspond,  respectively,  to 
the  propagators  for  the  empty  site  and  a  singly  occupied  site,  obey  a  set  of  coupled 
integral  equations.  Numerical  solution  of  these  equations  has  been  very  useful  in  ob¬ 
taining  quantitative  results  for  the  equilibrium  system,**  including  the  occupations 
of  the  two  spin-states  in  the  presence  of  a  magnetic  field.**  In  this  work  we  have 


generalized  the  non-crossing  approximation  to  non-equilibrium  to  produce  densities 
of  states,  occupations,  and  the  nonlinear  current  (2).  However,  as  a  large  spin- 
degeneracy  (large  N)  technique,  the  non-crossing  approximation  produces  a  Kondo 
peak  even  for  the  non-interacting  system  {N  =  1).  Consequently,  for  iV  =  2  in  a 
magnetic  held,  it  give  rise  to  spurious  peaks  in  the  density  of  states  at  the  chemical 
potentials. Therefore,  an  equations-of-motion  method  was  employed  to  comple¬ 
ment  the  non-crossing  approximation  and  isolate  its  shortcomings.  This  method 
corresponds  to  a  resummation  of  low-order  hopping  processes  emd  caimot  produce  a 
quantitative  description  of  the  Kondo  effect.  Nevertheless,  this  method  is  known^® 
to  give  the  right  qualitative  behavior  at  low  temperatures.  More  importantly  in  the 
present  context,  since  the  equations-of-motion  method  is  exact  for  iV  =  1,  it  gives 
rise  only  to  the  proper  Kondo  peaks  (as  identified  by  perturbation  theory**). 

The  equations-of-motion  method  consists  of  differentiating  the  Green  function 
G'{t)  with  respect  to  time,  thereby  generating  higher-order  Green  functions  which 
eventually  have  to  be  approximated  to  close  the  equation  for  G^{t).  The  procedure 
we  employ  here  is  the  same  as  the  one  used  in  Ref.  8,  which  in  the  infinite-f/  limit 


leads  to 


_ 1  -  M _ 

UJ  —  ~  ~ 
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with 


So<r(w)  =  2- 


and 


-  Cfcir  +  IT? 


/  \  _  fL/Rj^kf) _ 

"  ^■fc^K‘^-e.  +  e,-efc,  +  tV2r,’ 


(4) 


(5) 


where  /*/*(€)  ia  the  Fermi  distribution  in  the  left/right  lead  and  t,  is  the 
intermediate-state  lifetime.  g;{uj)  has  an  overall  amplitude  proportional  to  1  —  (n»), 
where  (n«)  is  the  occupation  of  the  other  spin-state.  Quantitative  calculation  of  the 
occupations  is  beyond  the  scope  of  the  equations  of  motion  in  the  present  approxi¬ 
mation  scheme.  Accordingly,  we  use  the  occupations  resulting  from  the  non-crossing 
approximation,  which  are  known  to  be  quantitatively  reliable  in  equilibrium.’^ 
Within  the  equations-of-motion  scheme,  the  Kondo  peak  for  spin  a  results  from 
the  self-energy,  Ej,(ci;),  due  to  virtual  intermediate  states  in  which  the  site  is  occu¬ 
pied  by  an  electron  of  opposite  spin,  a.  The  remaining  self-energy,  Eo^(u;),  is  the 
exact  self-energy  for  the  non-interacting  case.  Because  of  the  sharp  Fermi  surfaces 
at  low  temperature.  Re  {Eia.(u;)}  grows  logarithmically  at  u/  =  /it,*  ±  Ae,  giving  rise 
to  peaks  in  the  density  of  states,  —  4lmGj(u;),  near  those  energies.  The  peaks  for 
the  high-lying  spin  (low-lying  spin)  appear  near  u;  =  Pl.r  +  Ae  (w  =  pi,.n  —  Ac). 
At  zero-field  and  zero-temperature  equilibrium,  the  intermediate  states  giving  rise 


to  Ex^(u;)  have  an  infinite  lifetime,  and  the  true  peak  in  the  density  of  states  has 
an  amplitude  corresponding  to  the  unitarity  limit. Once  either  a  voltage  bias  or  a 
magnetic  field  is  applied  these  intermediate  states  acquire  a  finite  lifetime,  r^,  which 
cuts  off  the  logarithmic  divergence  of  Re  resulting  in  a  suppression  of  the 

peak  amplitudes.  The  lifetime,  r^,  of  spin  a  can  be  calculated  using  second-order 
perturbation  theory,  yielding 


a‘ 

X  [r^(e)rS(e  -  -h  e,0(l  -  A(e))/fl(e  -  e,  +  •  ( 


For  a  deep  level  at  zero  temperature  and  for  constant  T  this  simplifies  to 


_1_ 


1 

2irfi 


E 


rf,  0(mB  -  fiA  +  ^<r  -  e<r') 


fiB  -  fiA  +  ^  -  f*' 
{HA  -  -  C»')’ 


9‘' 


(7) 


which  explicitly  shows  that  the  lifetime  is  non-zero  only  for  finite  bias  or  finite 
magnetic  field. 

In  Fig.  1,  we  plot  the  density  of  states  for  two  spins  symmetrically  coupled  to 
two  leads,  consisting  of  Lorentzian  bands  of  width  2W,  so  that  r^((»;)  =  r^(w)  — 
rw^/2{u>^  +  W*),  with  r  =  1  and  W  =  100.  Results  are  shown  for  the  non-crossing 
approximation  (dashed  lines),  which  is  reliable  for  zero  magnetic  field,  and  for  the 


equations-of-motion  method  (continuous  lines),  which  has  the  correct  Kondo  peak 
energies  for  all  magnetic  fields.  In  equilibrium  and  zero  magnetic  field,  the  density 
of  states  exhibits  a  single  peak  at  the  Fermi  level  as  expected^^  (Fig.  la).  As  the 
chemical  potentials  split,  the  Kondo  peak  also  splits,  giving  rise  to  a  suppressed 
Kondo  peak  at  earh  chemical  potential  (Fig.  lb).  Upon  the  application  of  a  mag¬ 
netic  field,  the  densities  of  states  for  the  two  spins  become  different  and  the  Kondo 
peaks  shift  away  from  the  chemical  potentials  by  the  Zeeman  splitting  (Ac  =  0.2  in 
Figs.  Ic  and  Id).  The  peaks  move  up  in  energy  for  the  high-lying  spin  (Fig.  Ic) 
and  down  in  energy  for  the  low-lying  one  (Fig.  Id). 

The  main  conclusion  of  Fig.  1  is  the  emergence  of  new  energy  scales,  not  present 
in  equilibrium.  The  Kondo  peak  in  the  equilibrium  density  of  states  splits  out 
of  equilibrium  to  two  peaks  spaced  by  the  chemical  potentied  difference  A/t,  and 
suppressed  from  equilibrium  by  the  finite  dissipative  lifetime  In  Fig.  1(b),  the 
lifetime  broadening,  is  about  the  same  as  the  temperature.  It  is  apparent, 

however,  from  Figs.  1  (c)  and  (d),  that  neither  the  non-crossing  approximation  nor 
the  equations  of  motion  quantitatively  determine  the  Kondo  peaks  at  finite  magnetic 
field.  For  this  case  we  use  the  equations-of-motion  result  since  it  provides  a  good 
estimate  of  the  Kondo  peak  positions  (by  comparison  with  perturbation  theory). 
To  understand  the  shift  of  the  Kondo  peaks  with  magnetic  field,  it  is  helpful  to 


recall  how  the  peaks  in  the  density  of  states  derive  from  the  eigenstates  of  the 
system.  At  T  =  0,  G^{t)  involves  transitions  from  the  IV-particle  ground  state  to 
all  possible  iV  +  1  or  iV  —  1  states.  At  fl  =  0  the  correlated  ground  state  has  a  finite 
amplitude  to  have  an  empty  site,  and  thus  {c^.)  can  generate  trauisitions  from  the 
IV-particle  ground  state  to  the  ground  state  with  one  more  (one  less)  electron.  Since, 
by  definition,  the  ground  state  energies  differ  by  the  chemical  potential,  the  density 
of  states  includes  a  Kondo  peak  at  the  chemical  potential.  Within  a  variational 
calculation,^^  we  find  that  at  finite  magnetic  field  the  ground  state  is  polarized,  and 
adding  or  removing  an  electron  produces  no  overlap  with  the  new  ground  state. 
However,  there  is  a  correlated  excited  state  of  opposite  polarization  which  can  be 
reached,  and  which  consequently  gives  rise  to  a  peak  in  the  density  of  states,  shifted 
by  the  difference  in  energy  between  polarization  states,  i.e.  the  Zeeman  energy. 

The  current  follows  immediately  from  the  densities  of  states  (2).  In  particular,  the 
zero-temperature  current  is  the  integrated  density  of  states  between  the  two  chemi¬ 
cal  potentials,  weighted  by  the  coupling  to  the  leads  r,(u;).  At  zero  magnetic  field, 
therefore,  the  Kondo  peak  at  the  Fermi  energy  gives  rise  to  a  linear-response  conduc¬ 
tance  of  2e*//i  for  symmetric  barriers,  corresponding  to  perfect  resonant  transmission 
through  the  quantum  dot.“’^^  As  the  bias  is  increased  the  differential  conductance 
falls  rapidly  (Fig.  2a). This  occurs  firstly  because  the  differential  conductance  due 


10 


to  a  peak  in  the  density  of  states  must  fall  off  once  Afi  exceeds  the  peak  width,  and 
secondly  because  the  decreasing  dissipative  lifetime  suppresses  the  peak  amplitudes. 
Since  the  peadcs  in  the  density  of  states  persist  until  the  temperature  is  roughly 
one- tenth  the  coupling  to  the  leads,  F,  the  peak  in  the  differential  conductance  is 
observable  well  above  the  Kondo  temperature,  (Fig.  2a,  continuous  line). 

In  a  finite  magnetic  field  the  Kondo  peaks  are  shifted  away  from  the  chemical 
potential  so  they  contribute  very  little  to  the  conductance  in  linear  response.  As 
the  bias  is  increased,  however,  the  current  carrying  region  between  the  chemical 
potentials  grows,  until  at  =  Ac,  it  reaches  one  Kondo  peak  in  the  density  of 
states  of  each  spin  (see  inset  of  Fig.  2b),  In  Fig.  2b,  one  therefore  secs  peaks  in 
the  differential  conductance  at  A/x  =  Ae  (continuous  line).  In  fact,  by  comparison 
with  the  non-crossing  approximation  (Fig.  1),  we  expect  th";  equations  of  motion 
to  underestimate  the  full  strength  of  these  peaks.  Experimentally,  observation  of 
peaks  in  the  differential  conductance  at  A^  =  Ae  would  provide  a  “smoking  gun 
for  the  presence  of  Kondo  pnysics  in  transport  through  a  quantum  dot. 

In  this  work,  we  addressed  the  non-equilibrium  behavior  of  Anderson’s  model  for  a 
magnetic  impurity.  Experimentally,  the  model  describes  low-temperature  transport 
through  a  quantum  dot,  where  non-equilibrium  is  readily  accessible.  We  have  shown 
that  new  energy  scales  emerge  in  non-equilibrium.  Specifically,  the  difference  in 


n 


chemical  potentieds  and  the  inverse  dissipation  time  hlr^  lead,  respectively,  to 
splitting  and  suppression  of  the  Kondo  resonances  in  the  density  of  states.  Our 
results  have  led  to  a  novel  experimental  prediction  —  when  the  Zeeman  splitting  of 
the  spins,  Ac,  equals  the  applied  bias,  A^,  there  will  be  a  peak  in  the  differential 
conductance,  provided  these  energies  are  smaller  them  the  coupling  to  the  leads,  F, 
and  smaller  than  the  depth  of  the  levels,  Importantly,  this  signature  of 

the  Kondo  effect  persists,  for  a  wide  range  of  parameters,  to  temperatures  ~  r/20, 
which  may  be  magnitudes  larger  than  the  Kondo  temperature.  We  hope  that  this 
work  will  encourage  further  efforts,  both  experimental  and  theoretical,  to  probe  the 
non-equilibrium  physics  of  interacting  quantum  systems. 
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Figure  Captions: 

(1)  Density  of  states  for  an  Anderson  impurity  symmetrically  coupled  to  two  leads 

with  chemical  potentials  and  /Xfl(=  0)  and  Lorentzian  bandwidth  2W,  from  the 
equations-of- motion  method  (continuous  line)  and  the  non-crossing  approximation 
(dashed  line).  The  impurity  has  two  spin  states  with  energies  Cf  and  and  an 
on-site  interaction  U  —*  oo.  All  energies  aire  in  units  of  the  total  coupling  to  the 
leads,  r.  The  band  width  is  W  =  100  and  the  temperature  is  T  =  0.005,  roughly  a 
factor  of  two  lower  than  the  magnetization  Kondo  temperature  [Ref.  14).  (a)  The 
equilibrium  =  0)  density  of  states  at  zero  magnetic  field  =  —2.0.  The 

density  of  states  exhibits  a  single  peak  at  the  Fermi  level  [Ref.  17).  (b)  The  non¬ 
equilibrium  {fic  =  0.3)  density  of  states  at  zero  magnetic  field  e,  =  =  —2.0.  There 

is  a  suppressed  Kondo  peak  at  each  chemical  potential.  (c),(d)  The  non-equilibrium 
(^^  =  0.3)  density  of  states  for  spin  up  (c)  and  spin  down  (d)  at  fimte  magnetic  field 

=  -1.9,  Cl  =  -2.1.  The  Kondo  peaks  shift  away  from  the  chemical  potentials  by 
the  Zeeman  splitting  Ac  =  0.2;  the  shift  is  up  in  energy  for  the  up  spin  and  down 
in  energy  for  the  down  spin. 

(2)  Differential  conductance,  e  dJ/dA/x,  with  hr  =  0,  vs.  applied  bias,  (a)  Zero 

magnetic  field  differential  conductance  via  the  non-crossing  approximation,  (b)  Dif¬ 
ferential  conductance  at  the  finite  magnetic  field,  Ac  =  0.2,  used  in  Figs.  1  (c)  and 


i.l 


(d),  via  equations  of  motion.  As  shown  in  the  inset,  when  the  chemical  potential 
difference,  A/i,  reaches  the  Zeeman  splitting,  Ac,  the  Kondo  peaks  in  the  density  of 
states  enter  the  region  between  the  chemical  potentials,  giving  rise  to  a  peak  in  the 
differential  conductance. 
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